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Abstract 

Nonperturbative infrared finite solutions for the gluon polarization ten- 
sor have been found, and the possibility that gluons may have a dynamically 
generated mass is supported by recent Monte Carlo simulation on the lattice. 
These solutions differ among themselves, due to different approximations per- 
formed when solving the Schwinger-Dyson equations for the gluon polarization 
tensor. Only approximations that minimize energy are meaningful, and, ac- 
cording to this, we compute an effective potential for composite operators as 
a function of these solutions in order to distinguish which one is selected by 
the vacuum. 
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During the last years there has been much effort in trying to obtain a non- 
perturbative form for the gluon propagator [1-8] , and perhaps one of the most inter- 
esting results is the one where it is argued that the gluon may have a dynamically 
generated mass [[|. The existence of a mass scale, or the absence of a pole at k 2 = 
is natural if one assumes that gluons do not propagate to infinity, i.e. these prop- 
agators describe confined gluons. From the phenomenological point of view, this 
possibility may shed light on several reactions where long distance QCD effects can 
interfere, and examples about the consequences of massive gluons can be found in 
the literature, see, for instance, Ref. [9-11]. 

The study of the infrared behavior of the gluon propagator was also performed 
numerically fll2"l , and the last numerical lattice computations give strong evidence 
for an infrared finite gluon propagator [13]. The most recent results |13[ seems to 
exclude a propagator of the form 1/k determined in several analysis of Schwinger- 
Dyson equations (SDE) and claimed to give confinement Jy, £| 3- In particular, 



the results of Bernard et al. [[T3|] are consistent with an infrared finite propagator, 
but can be fitted by the Cornwall's massive propagator p3] as well as by the one 
found by Stingl et al. |5J. These lattice computations could be criticized because of 
the gluon propagator gauge dependence. However, there is also a lattice study of 
the "gauge-invariant" two-point correlation function of the gauge field strengths by 



Di Giacomo et al. []14 |, where a finite correlation length for the gluons in the QCD 
vacuum is found, indicating the presence of a mass scale for the gluon propagation. 

If the lattice results can be considered exciting, the same cannot be said about 
the analytic calculations. Looking at Ref. [l[] to || we verify that infinite as well 
as infrared finite non-perturbative propagators have been found in the literature, 
although this difference can be credited to the fact that some authors discarded 
from the beginning an infrared finite solution for the gluon propagator. Even if we 
separate the solutions of SDE for the gluon self-energy, according to the existence 
or not of a pole at k 2 = 0, we still remain with different forms of non-perturbative 
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propagators, which reflect the different approximations involved in the calculations, 
for instance, different truncations of the integral equations, or the handling of higher- 
point functions which are present in the two-point one, etc... 

In this work we propose to compute the effective potential for composite op- 



erators |T5[ at stationary points as one method to distinguish among the different 
solutions, i.e. we will calculate the vacuum energy for some infrared finite non- 
perturbative Schwinger- Dyson solutions of the gluon propagator, and the vacuum 
will select the solution that leads to the deepest minimum of energy. Therefore, this 
will be a complementary tool to find solutions of the gluonic SDE, i.e. they should 
be tested by means of the effective potential calculation. As we are far from being 
able to solve the SDE without any simplification, only approximations that mini- 
mize the vacuum energy should be performed when solving these integral equations. 
We will also discuss how the SDE solutions can be constrained phenomenologically. 
For a non-abelian gauge theory the effective potential has the form 



V(S,D,G) = -ijJ^TrQnS^S-^S + l) 

~ % I J^ Tr{lnG ° lG " G ° lG + 1} 

4/(& Tr(ln ^"^ + 1) 
+V 2 (S,D,G), (1) 

where S, D and G are the complete propagators of respectively fermions, gauge 
bosons and Faddeev- Popov ghosts; So, D Q and Go the respective bare propagators. 
V 2 (S, D,G) is the sum of all two-particle irreducible vacuum diagrams, depicted in 
Fig.(l), and the solutions of the equations 

^ = ^ = ^ = (2) 
SS 5D 5G ' K ; 

give the SDE for fermions, gauge bosons and ghosts. 
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We can represent 1^(5", D, G) analytically by 

iV 2 (S,D,G) = -~Tr{TSTSD)- -Tr{FGFGD) 

+hrr(V^DV^DD) + ~Tr(r^DD) (3) 

where T, F, and are respectively the proper vertex of fermions, ghosts, 
trilinear and quartic gauge boson couplings JTB|. In Eq.( |3|) we have not written the 
gauge and Lorentz indices, as well as the momentum integrals. 

The complete gauge boson propagator D is related to the free propagator by 

zr 1 = Dq 1 - n, (4) 

where II is the gluon polarization tensor, which is obtained from Eq.([l]) and Eq.(||]), 
and described by 

II = TSTS + FGFG - -T {3) DT {3) D - -T W D. (5) 

2 2 K ' 

The diagrams contributing to II are shown in Fig. (2), and this self-energy is the one 
that has been solved in many different ways, leading to a series of non-perturbative 
forms for the gluon propagator in the infrared region [1-6]. As our intention is to 
compute the effective potential through a variational approach, the best ansatz for 
the complete propagator will be given by the approximate solution of II. 

The physically meaningful quantity we need to compute is the vacuum energy 
density given by 

Q = V(S,D,G)-V(S P ,D P ,G P ), (6) 

where D p is the perturbative counterpart of D, and we are subtracting from the po- 
tential the perturbative part that does not contribute to dynamical mass generation, 
denoted by V(S P , D p , G p ). As our interest is concentrated in the pure glue theory, 
we disregard all the quark contributions for the effective potential from now on, i.e. 
we drop the first contribution on the right-hand side of Eq.( |3[). It is opportune to 
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recall that for perturbative gluons is plagued by infrared divergences, and only for 
the infrared finite gluon propagators it can be calculated without ambiguities ||. 

Before we start calculating Q it is convenient to discuss how the solutions of 
II are determined. One of the first calculations of the gluon polarization tensor is 
due to Mandelstam [I] , where the ghosts and the diagram with quartic coupling are 
disregarded. The neglect of ghosts diagrams was shown to be correct, because their 
contribution is numerically small PJ. These approximations, with the use of the 
Landau gauge, were shown to be satisfactory in the lengthy and detailed work of 
Brown and Pennington ||J]. With this approximation only the diagram with trilinear 
coupling should be considered in Eq.( |3|), i.e. the gluon polarization tensor is going 
to be given by 

n = At^dt^d. (7) 

More complex approximations can be used, but in practice only the diagram with 
the trilinear coupling has been considered in most of the works up to now. As far 
as we know, few papers have gone beyond these approximations, and one of these 
approaches is the Cornwall's determination, through the use of the pinch-technique, 
of a set of diagrams leading to a gauge invariant SDE ||. Another one is the Stingl 
et al. manipulation of the higher-point functions |J, although this calculation was 
done in Landau gauge. 

The vacuum energy can be computed with the same approximations performed to 
obtain solutions of the gluon self-energy. However, a much better approximation to 
study the vacuum stability results when we plug the stationary conditions (Eq.( |^) 
and Eq.(^|)) into Q, obtaining an expression for the values of Q at its minimum, 
which will be denoted by (Q). For example, with the Mandelstam approximation 
the 2PI diagrams contributing to fl would be given by 

tt 2 = —Tr(T (3) DT {3) DD) - (D -> D p ), (8) 
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but applying the stationary condition Eq.(^) we obtain 

(n 2 )= t -Tr(UD)-(D^D p ), (9) 

where the calculation is now reduced to an "one-loop" integration. It is important to 
stress that (O) is much less dependent on the ansatze we use to compute the vacuum 
energy. This has already been emphasized in calculations of the vacuum energy for 
chiral symmetry breaking |T7|]. Note that (O) also does not suffer from possible 



problems in the determination of the effective potential for composite operators [|T^] . 

The vacuum energy at stationary points, within the approximation discussed 
above, is given by 

(O) = t -Tr[\n(D l D)-D 1 D + l} 

-\tt[{D^ - D^){D - D )l (10) 

where we assumed D p = D . 

Displaying the momentum integration, using the Landau gauge expression for 
the complete gluon propagator 

d^(p 2 ) = - 2 * 8 . or - ^f), (ii) 

p z — LL{p z ) p z 
and calculating the trace, Eq. ([T0|) can be cast in the following form 



= _3(jV 2 -D f *P 



(27T) 



n , n v 2 n 2 

ln(l + — ) + 



(12) 



p 2 + n v p 2J 3P 2 (P 2 + n) 

where all the quantities are in Euclidean space (P 2 = — p 2 ), and N = 3 for QCD. 
Note that the computation of the full vacuum energy Q requires the use of nonpertur- 
bative vertices satisfying the Slavnov- Taylor identities, which give strong relations 
between the diagrams of Fig.(l) ||19|| , and we escaped from this difficulty in our (Q) 
calculation. 

To compute (Q) we need only a reasonable ansatz for the gluon propagator 
or its self-energy EE, which will be given by the SDE solutions to be presented in 
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the sequence. One of the infrared finite propagators found in the literature was 
determined by Cornwall || 



D-\P 2 ) = [P 2 + m 2 g (P 2 )}bg 2 ln[ + A / 1 '- ], (13) 
where m g (P 2 ) is the dynamical gluon mass 



m 2 Q (P 2 )=m 2 g 



( P 2 + 4m 2 , \ Am 2 



-12/11 



(14) 



a 2 /' a 2 

m s is the gluon mass scale, g s the strong coupling constant, b = 33/487T 2 is the 
leading order coefficient of the f3 function of the renormalization group equation, 
and A = Aqcd is the QCD scale. This form for the propagator was obtained as a fit 
to the numerical solution of a gauge invariant set of diagrams for the gluonic SDE, 
which also imposes g s ~ 1.5 — 2. 

Eq.([L3|) is consistent with the expected ultraviolet behavior of the gluon propa- 
gator 

^(P 2 )^P 2 ln(P 2 ). (15) 

However, the dynamical gluon mass m g (P 2 ) is expected to decrease faster asymp- 
totically as predicted by the operator product expansion (OPE) pO], ||. It should be 
remembered that the SDE for the gluon propagator is accurate in the determination 
of the infrared dynamical mass, but its ultraviolet behavior can be better determined 
by means of the OPE. The behavior of the dynamical gluon mass through the OPE 
was determined by Lavelle ||, and is equal to 

TT (P 2 \ 34jV7r (a s G 2 ) 

where (a s G 2 ) is the gluon condensate. Note that as long as the gluon self-energy is 
infrared finite, and no matter its ultraviolet behavior is given by Eq.flU]) or Eq. (|T6|) , 
it is easy to see that (Q) is free of divergences 0. 

Another infrared finite propagator has been found by Stingl et al. 0, its form 
agrees with that derived by Zwanziger based on considerations related to the Gribov 
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horizon |7]], and is given by 

B -< p2 > = • < 17 > 

r (P 2 +c 2 ) 

where fi s is a mass scale not determined in Ref. 0. Eq.( |T7p with c = is the one 
found by Zwanziger |jj . As claimed in the last paper of Ref. || , the full propagator 
should contain a c ^ term, however, due to the complexity introduced by such 
factor, this case was not studied in detail. The term \x\jP 2 is exactly what is 
expected by OPE analysis whenever a mass scale for the gluon is introduced, as 
shown by Eq.( [16]), therefore, its appearance is not surprising. With c = the 
propagator goes to zero at P 2 = 0, however the gluon self-energy diverges, and we 
discard this possibility here assuming that c = fi s , which could be expected as a 
natural mass scale leading to a finite gluon polarization tensor. The Bernard et 
al. lattice result for the gluon propagator (T3| can be fitted by Eq.flTSP as well as 
Eq.flU- 

Finally, Marenzoni et al. also performed a lattice study of the gluon propa- 
gator in the Landau gauge, obtaining for its infrared behavior the following fit 

D ^ K ") = 2 ( 18 ) 

mf + ZP 2 {jj)v 

where mi, Z and 77 are constants determined with the numerical simulation, mi is of 
0(A ~ 160 MeV), Z ~ 0.4, and rj ~ 0.5 what is slightly different from the previous 
propagators. The results of Bernard et al. also show the behavior (-P 2 ) 77 , but with 
a smaller value for rj. We will consider only the propagators discussed up to now. 
Although other propagators can still be found in the literature, the propagators that 
we presented above are the most discussed ones and enough for our purposes. 

The result of Ref. 0, which led to Eq. (pT3[) , was shown to be gauge invariant, 
and the others were obtained in the Landau gauge, which we assume here expecting 
that any gauge dependence should be negligible. In the above equations for the 
infrared propagator, except from Cornwall's one, it is not clear how to connect these 



expressions with their ultraviolet part, and this is an extra complication, because 
any, ad hoc, link between the ultraviolet and infrared propagator will introduce a 
new mass scale. It is obvious that the physics should not depend on this particular 
scale, but it does introduce some arbitrariness. We, therefore, will not consider 
the ultraviolet renormalization group logarithmic behavior of the full propagator, as 
described by Eq.([L^), since it should be common to all propagators (although it is 
only present in Eq.flTB])). In the especific case of the lattice propagator, Eq.([L8|), in 
which there is no clear separation between the infrared and ultraviolet regions, we 
shall simply neglect its peculiar form, assuming that the vacuum gluon polarization 
can be cast as a pole in the perturbative propagator factorized from Eq. (|T8D . 

We adopt the following ansatze to calculate (Q): a) IT(P 2 ) = m 2 (P 2 ), where 
m 2 (P 2 ) is given by Eq. (|14D ; b) To be consistent with the massive Cornwall propa- 
gator and with OPE, we will consider a gluon self-energy that interpolates between 
the constant infrared behavior of Eq.(|l4|) and the ultraviolet one of Eq . (|l6t) , which 
is given by 

n(p 2 ) = ^ 2 -^ 2 ) + ^(^-/i 2 ); (19) 

c) II(P 2 ) consistent with Eq.([T7|), i.e. 

n(P 2 ) = ^ (20) 

d) We assume that IT(P 2 ) for the lattice inspired propagator is given by 

II(P 2 ) = - P 2 ) + j^fj^OiP 2 - (21) 

which follows from considering the generated gluon mass as a pole in the perturbative 
propagator, factorized from Eq . fllS|) . Finally, since m g and fi s are not determined 
by the solution of the SDE, we compute the vacuum energy as a function of these 
parameters. 

In Fig. (3) we show the vacuum energy computed with our ansatze a) to d), as 
a function of the gluon mass scale (m g = fi s = mi). Note that once selected a 
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form for the gluon self-energy, II, the (Q) calculations showed in Fig. (0) only give 
us the vacuum energy as a function of the gluon mass scale (or m g /A), and we do 
not have a criterion to fix a particular value of this scale parameter. We can only 
compare the different vacuum energy curves, corresponding to the different gluon 
self-energy behavior, assuming a common value for the gluon mass scale. From the 
behavior of the vaccum energy presented in this figure we verify that the Cornwall 
propagator, Eq. fllTf) , leads to the deepest values of minimum of the vaccum energy, 
as long as m g > 1.6 A, compelling us to believe that this could be the most favoured 
form of the gluon propagator to match the one chosen by the QCD vacuum, for 
gluon masses obeying such relation. For smaller gluon masses {m g < 1.6A), the 
lattice propagator is the one that best represents the actual gluon propagator. It is 



clear that apart from the lattice propagator (Eq.(18)), where all the constants are 
determined, we can only compare the different values for the vacuum energy, if we 
are able to constrain the parameters m g ,[i s . In principle they should be related to 
the scale A, but we are far from this achievement, and we discuss in the sequence 
some of the constraints that have been put forward on these parameters. 

The phenomenological constraints on the dynamical gluon mass may come from 



the quarkonium decays into gluons |§ |TT| , potential models for glueballs [pi] , ap- 
proximations to the QCD-Hamiltonian combined with a variational method [p^| , a 
possible relation between the gluon condensate and the gluon mass through QCD 



sum-rules ||23|| , the two-gluon exchange contribution to the Pomeron ||10|| , and the 
vacuum functional of SU(3) Yang-Mills theory ||24|| . The gluon masses obtained in 
most of these studies are of the order of 2A, but one should be careful when consid- 
ering this value because with few exceptions these calculations consider only a bare 
gluon mass, and the specific behavior of the dynamical mass with the momentum 
may modify the result. These calculations in general also involve many different 
(and some untenable) approximations. Therefore, we are in the safe side if we con- 
sider that the gluon mass may lay between A and 3A, and according to Fig. (3), the 
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propagator determined by Cornwall and the one obtained in the lattice simulation 
lead to the deepest vacuum energy, and approach the actual behavior of the con- 
fined gluon propagator. Finally, we can also use one relation found by Cornwall [|3| 
between the gluon condensate and the vacuum energy 



(g*G%G?) = < fi > > ( 22 ) 

which allows a determination of the gluon mass. To arrive at this relation it has 
been assumed that the f3 function is given by its perturbative behavior (/3 ~ —bg 3 ). 
Although the f3 function can be different from this in the condensed phase, maybe 
this relation is the most trustable one to allow a determination of the gluon mass. 



The left-hand side of Eq. ( p!q) has been evaluated through QCD sum-rules [f25| and 
is equal to 0.47 GeV 4 . Therefore, assuming A = 0.3 GeV we obtain the following 
gluon masses in the cases we discussed above: a) 527 MeV, b) 627 MeV, c) 715 MeV, 
and d) 560 MeV. Considering all the gluon mass determinations in the literature we 
verify again that the propagators of Ref. |J and Ref. |TB[ are the ones to provide 
the most consistent picture. 

In conclusion, starting from the effective action for composite operators, we 
determined an expression for the vacuum energy of QCD without fermions in the 
Landau gauge. The vacuum energy calculation was used as a criterion to determine 
which, among a series of different forms determined for the gluon propagator, leads 
to the deepest minimum of energy. We verify that two expressions are clearly 
favored and, consistently, when we compare the vacuum energy with the gluon 
condensate, we obtain gluon masses for these propagators that are in the range of 
several determinations of this mass existent in the literature. 
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V 2 (S,D,G) = 




Figure 1: Two-particle irreducible vacuum diagrams contributing to the effective 
potential. 





Figure 3: Result of the vacuum energy calculation as a function of the gluon mass 
for different ansatze for the gluon vacuum polarization. The curves are related to 
the ansatze assumed inside the text: a) solid line Eq. (|14D ; b) dotted line Eq. (|19D ; c) 
dashed-dotted line Eq.@g); d) dashed line Eq.(fZJ). 
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